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Abstract 

Here we address the problem of bosonizing massive fermions without making 
expansions in the fermion masses in both massive QED2 and QED% with iV fermion 
flavors including also a Thirring coupling. We start from two point correlators 
involving the £7(1) fermionic current and the gauge field. From the tensor structure 
of those correlators we prove that the ?7(1) current must be identically conserved 
(topological) in the corresponding bosonized theory both in D = 2 and D = 3 
dimensions. We find an effective generating functional in terms of bosonic fields 
which reproduces those two point correlators and from that we obtain a map of the 
Lagrangian density tp r — m)ip r into a bosonic one in both dimensions. This 
map is nonlocal but it is independent of the eletromagnetic and Thirring couplings, 
at least in the quadratic approximation for the fermionic determinant. 

PACS-No.: ll.15.Bt , 11.15.-q 
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1 Introduction 



One among many dreams of the theoretical physicists nowadays is the possibility of ex- 
tending to higher dimensions (D > 2) the bosonization of fermionic models. This can 
be justified by some good properties of one formulation in contrast with the other one. 
For instance, strong coupling physics in one model corresponds to weak coupling in the 
other one. A classic example is the map between the massive Thirring and Sine-Gordon 
models in D = 2. Another interesting aspect of such map is the fact that the usual 
electromagnetic charge in the fermionic model corresponds to the topological charge of 
the associated soliton field. Furthermore, we can have a map between a linear theory like 
massive free fermions and a nonlinear one (Sine-Gordon at f3 2 = Air) on the other side. 

In view of these and other interesting properties a lot of work has been devoted to 
the issue of bosonization 13 El ( see also [I] ). There has also been many attempts to 
generalize those ideas to higher dimensions [5-22]. For massive fermions in D = 2 most 
of the methods is based on expansions around massless fields which are local conformal 
theories. In D = 3, although the case of massless free fermions can still be mapped 
into a bosonic theory ( see [H] ) such theory is nonlocal. Besides, the conformal group is 
finite in D = 3 and not so powerful as in D = 2 which makes expansions around the 
massless case nontrivial. The other possibility is to employ functional methods. Once 
again the case of massless fermions is easier to deal with since the fermionic determinant 
can be exactly calculated for m = 0. For massive fermions in D = 2 a nontrivial Jacobian 
under chiral transformations plays a key role in deriving the Sine-Gordon model ( see [26- 
30]). In D = 3 , chiral transformations play no role and although some nonperturbative 
information is known [21] about the fermionic determinant we are basically left with 
approximate methods like the one used in this work. On one hand, the method used here 
is inspired in the approach carried out in [19] which is rather simple and based on two 
point correlation functions. On the other hand, our bosonization rules departure from 
those in jT^] in the sense that they are independent of the interactions. We extend the 
approach of 19] by introducing an eletromagnetic interaction and making use of the 1/N 
expansion which allows us to go beyond the lowest-order in the coupling constants. In the 
next section we introduce the model we are working with and obtain a general expression 
for the generating funtional of the current and gauge field correlators. The expression is 
valid for arbitrary dimensions and depends on the vacuum polarization tensor. In sections 
3 and 4 we make the calculations explicit in D = 2 and D = 3 dimensions respectively. We 
first obtain in the fermionic theory the two point current correlation functions involving 
also the gauge field and then we write the current in terms of bosonic fields and derive 
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the corresponding action for such fields which reproduces their correlators. In the final 
section we draw some conclusions and comment on similar approaches in the literature. 



2 Generating Functionals 

We start by introducing the notation which will be used in both D = 2 and D — 3. The 
generating functional for a generalized QED with Thirring self-interaction is given by 



V K 



N 

= /^,n 

J r=l 



Vtp r Vipr exp lij d D x 



4)A 



^ ( Vvt^Vv) 2 + \ {d^f + j; (^A) + k, a* | , (l) 



where N is the number of fermion flavors and summation over the repeated flavor index 
r (r = 1,2, ...,N) is assumed. It is convenient to introduce an auxiliary vector field 
and work with the physically equivalent generating functional: 



j;, K„] = J VApVBp V^j r Vi> r exp^i J d D x 



(2) 



'N VN 
Integrating over the fermionic fields we obtain 
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(3) 



Now we compute the fermionic determinant keeping terms of order (1/N)° and 
(1/iV) 1 / 2 . Higher order terms will be neglected. Furthermore, since we are only inter- 
ested in two point correlators, the terms higher than quadratic in the sources will 
not be taken into account either. This amounts to the quadratic approximation for the 
fermionic determinant which gives 

j;,K,} = JvA.VB.exp'-J £^ {-A,[e^k 2 {I - X) + \k 2 g^] A u + B»B„ 



N A R A R 
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where d^ v = g^ u — ^Jf- and the tilde over the fields represent their Fourier transformations 
in momentum space. The quantity IP" is the vacuum polarization tensor: 

n ""<*>= *■/(§> 

In order to proceed further we have to calculate IP" which depends on the dimension- 
ality of the space-time. 
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•ft — m + ie (ft + tf.) — m + ie 



(5) 



3 Bosonization from two point correlators in D = 2 



In this section we restrict ourselves to the D = 2 case. Using dimensional regularization 

4m? 



we obtain, below the pair creation threshold (z = -J^p < 1), : 



with 
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; < z < 1, 



z < 0. 



(7) 



(8) 



2 \A^-i) vv^-^-v-*, 

Once the tensor IP" is calculated one is left with a gaussian integral over the vector fields 
and from which a generating functional quadratic in the sources is derived. Such 
generating functional furnishes the following two point correlators 1 



TT 2 (p 2 - k 2 n 2 ) 

<f»(k)j s A-k)> = ^ g V + nrr (9) 

<3l{k)A v {-k)> = <r, (10) 

\k 2 D ) Afc 2 

1 More precisely, we should have written explicitly the two point functions in the form < G(k)H(p) > = 
I(k)S^ (k+p) but in this article we will not display the delta function for a matter of convenience. Notice 
also that when we write A rs = F + G8 rs it is assumed that F multiplies a N x N matrix where all entries 
are equal to one. 



<A»(k)A"(-k)> = (-_ + ^_ r- + f-, (ii) 
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where 

D = [ne 2 -A; 2 (l + <7 2 n) 

The tensor structure of the above correlation functions are in full agreement with the cor- 
responding Ward identities based on the U (1) symmetry and it will play a key role in our 
bosonization procedure. At this point it is important to stress that our approach deviates 
from that of reference [HI]. The 1/N expansion we have relied upon, which coincides 
with the quadratic approximation for the fermionic determinant, is not equivalent to the 
second order weak coupling expansion used in |19j . 

Now, in order to derive a bosonized expression for the currents j r = ip r ^^ipr, we write 
down the most general decomposition for a vector in the momentum space: 



fp{k) = epsk s <j> r {k) + k^ r (k). (13) 
Substituting it in and using the identity e^sk & e ai k' 1 = k 2 6 a p we obtain 

< f a (AO j' p i-k) > = -k 2 e af3 < <\? (k) <p s {-k) > -k a k p < <p r (k) <p s {-k) > 



, 5 k s kf3 < r (k) tp s (-k) > -efrWhc < (p r (k) <p s {-k) > 



n 2 (e 2 - k 2 g 2 ) 

From the above it is not difficult to derive 

< ip r (k) if 3 (—k) > = 0, 

< <f) r (k) v s i-k) > = 0, 

< <P T (jfe) <p s {-k) > = 



@a/3 + n^Q.^^ 



a/3 < 



1 

k 2 



n 2 (e 2 - k 2 g 2 ) 



+ n 5 T 



(14) 

(15) 
(16) 

(17) 



from which one can safely set (p — 0. On the other hand, substituting a general decom- 
position 



< (AO A„ (-AO > = M e^s k s + Qk 
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in the mixed correlation functions < j r (k)A v (—k) > given in (fTU|) we conclude that 



<f(k)A^-k) > 



e n 



Nk 2 D 



e^sk 6 



(19) 



Now we are in a position to derive the bosonic Lagrangian density CsiA^, cff) which 
is compatible with the correlation functions (|ll )l .(|17 jl and p9j) . For this purpose we start 
from the following Ansatz 



C B (A 



I- 1 ' 
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(20) 
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where R rs , S r , T\ and T 2 will be determined as follows. We introduce the external sources 
X r and and define the generating functional 

r N 

Z B [X r ,K^ = / f[ V( pr VA ^JS X [C B (A^) + X^ + K^]^ (21) 
J r =l 

Assuming that R rs is a symmetric non-singular matrix we have performed the gaussian 
integrals in (}2*Tj) and obtained an explicit formula for Zb [X r , K^] from which the two- 
point correlators can be obtained. By matching these correlators with (|TT|) . (fT7j) and (fTT?|) 
we determine the bosonic Lagrangian uniquely : 



■ N k 2 2 I N \ / N 



A a A 13 



(1 - \)k 2 9 a p + \k 2 g a t 



/ N 
\r=l 



[22) 



All the steps which lead us to (j22j) are technically simple and not very elucidating. We 
should only mention that no expansions on either 1/N or any of the coupling constants 
have been made in those intermediate calculations. Furthermore, notice that if we quan- 
tize Cb{A, Ji , r ) and integrate over the scalar fields <p r in ()22|) this will lead to a nonlocal 
effective action for the photon which was studied in [21] where we concluded that, al- 
though nonlocal, the theory is free of tachyons. Next, by comparing with the Lagrangian 
density written in terms of fermionic fields in (1) we have the bosonization formulae for 
each fermion flavor (no sum over repeated roman indices below) : 



W7mVv(A0 = e^k v (f> r (k), (23) 

-&(*)(Jf + m)Vv(*) = \f^ r - ( 24 ) 

Now some comments are in order. First of all, if for some given flavor we do a U(l) 
transformation {ip r —>■ e ta ip r ) in the expectation value < j r > and use any regulariza- 
tion scheme preserving the U(l) symmetry it will be easy to derive the Ward identity 
< j r Jd v jl >= which implies the tensor structure < j^jl > oc 6^ u and consequently 
we will have (J23|) . So the current is topological due to the U(l) global symmetry and 
that must hold non-perturbatively. On the other hand, the bosonization rule ()24j) is only 
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approximate since the full expression would require, in our approach, the complete knowl- 
edge of the fermionic determinant which is only possible for m — > 0. In this case II — > 1/tt 
and we end up with N massless scalar fields topologically coupled to the gauge field. 
Integration over the gauge field leads to N — 1 massless scalar modes and one mode with 
m 2 = e 2 / (Nir + g 2 ) . Thus, reproducing the particular case of the so called Schwinger- 
Thirring model for N — 1 j23|, as well as the Schwinger model result m 2 = e 2 /7r for g — > 
and N — 1. In the opposite limit of large mass (z — > 0) we have tt — > —2z/3 at leading 
order. Substituting back in (}2*2*|) we arrive at a divergent result at m — > oo which is in 
agreement with the m — > oo limit of the corresponding Sine-Gordon model. See [22] for 
a similar comparison in the case of the massive Thirring model without eletromagnetic 
coupling. 

4 Bosonization from two-point correlators in D = 3 

In D = 3 dimensions the vacuum polarization tensor (0) calculated by means of dimen- 
sional regularization is given by : 

W u {k) = i Ut E" u + U 2 k 2 6^ , (25) 
with = e^ up k p and, in the range < z < 1, 




ili = -7- In —- : n 2 = 

87rzV2 \l-zV 2 )' 16nmz 

While for z < we have 
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(26) 



n! 



arctan • 
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Wnm z 



1 - 



l + z 



2(-z) 



1/2 



arctan 



(27) 



4vr(-^) 1 /2 

In fact the 111 amounts to a regularization dependent finite term [313, which was taken 
equals to zero due the dimensional regularization used. 

Substituting (|23|) in the general expression (J3J) we can obtain the two-point functions: 



<fMjl(-k)> 



<f(k) A u {-k)> 



1 

'N 



i, ( D ' + 3 ( " +iii 'K)^ 



e 



[P6^ - i^E, 



(28) 
(29) 



7 



< (k) A u (-k) > 



9 flu 
Xk 2 
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(30) 



where we found convenient to define 



P 
Q 



e 2 -k 2 g 2 )(k 2 Tl 2 .-Tli)-k 2 Tl 2 



e 2 - k 2 g 2 )U 2 - 1 



- e" 



(31) 
(32) 



In Analogy to the D = 2 case, we use for the U(l) current a general decomposition in the 
momentum space : 



f a (k) = eafrkP B?(k) + k a( p r (k), 



(33) 



from which we get 



<fa(k)fp(-k)> = -k a kp<<p r (k)<p s (-k)>+E ai E 5 p<B~J(k)B 5 s (-k)> 

+ k a Eps < r (k) B 5 S (-k) > +k/3 E ai < B? (k) <j) s (-k) > .(34) 

Multiplying the last expression by fc a Ar we conclude that < cff (k) <p s (—k) >= 0. Now, 
multiplying the resulting expression by k a we have E a p < <p r (k)B^(—k) >= 0. A similar 
manipulation was used in the last section to derive (JT3J) , (HHJ) and (fTTj) . Concluding, we 
can certainly neglect the scalar fields r = and minimally bosonize the U(l) current in 
D = 3 in terms of a bosonic vector field. The bosonic version of the current is once again 
of topological nature and identically conserved. As in D = 2 case this happens because of 
the U(l) global symmetry of the fermionic Lagrangian. Taking cff = and substituting 
the decomposition ()33)1 in (f2*KJ) and (|2l?|) . after some trivial manipulations, we end up with 

<B?(k)B s s (-k)> = -C rs g' 5 + (c - ff^ s -^^EP s , (35) 

<A'(*,B?(-*)> = -(^ + D ) r -(jy}E»» + D,^, (36) 

where C rs and D s are arbitrary and 



{H\)rs 
{H2)rs 



(37) 
(38) 
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The arbitrariness of C rs and D s shows that the bosonization fields must be gauge 
fields but the corresponding gauge symmetry is independent of the eletromagnetic one. 
Analogous to the last section, we next derive a Lagrangian density compatible with the 
two-point correlators . and (J3Uj) . Let us suppose we have a bosonic Lagrangian 
density £#(^4^, B") of the form 

C B (A^, B v r ) = B? O r ; u B v s + u'B» Ep, A u + A„A„ ( Vl 6^ + v 2 g^) , (39) 
Where 0™ has a general tensor structure in both momentum and flavor space, 

O r ;„ = oT s 0^ + b rs E^ + d rs 9fiV (40) 

The quantities a rs ,b rs ,d rs ,u s ,Vi,V2 will be determined by matching the two-point corre- 
lators derived from the generating functional : 




with the correlators (|3*3j) , ()36|) and (J3*U|) . Assuming that the matrices a rs , b rs , d rs are sym- 
metric and non-singular, our first step is to perform the gaussian integrals over the N 
vector fields B v r which furnishes: 



Zr \Y r K 1 = f DA e* I dH-K Y f+ u rApE^)(O-%t(Y»+u s E"0A,3)+A»A»(vi6^+v^ 

(42) 

The inverse operator 0~ x has the same structure of (jJUJ): (C _1 )^t = a rs @nv + b E^ v + 
<^ rs 5V wnere the tilde variables are functions of the non-tilde ones. In particular, d rs = 
(d _1 ) rs . Since we have the contractions E^ u E ua = —k 2 9£; E fiv ua = E£ it becomes clear 
from (}4*2"j) . with a little thought, that even after integrating over A^ the only term in log Zb 
quadratic in the sources YJ 1 and contracted via the metric tensor is : — \ Y^Y^d rs g PiV . 
Therefore the two point functions will require the identification C rs = (d/2) rs = 
(d-V2)". 

Given that the matrix C rs is arbitrary, in order to simplify the calculations, we have 
chosen it proportional to the identity and consequently we reduce the undetermined pa- 
rameters in our original Ansatz for the bosonic lagrangian since d rs = dS rs . For analogous 
reasons, after integrating over A^, the crossed terms in the sources YJ 1 and K v will be 
only contracted by tensors E^ v and 6^ u . Thus, the term proportional to g u ^ L on the right 
handed side of the correlation functions ()36|) must vanish, which fixes D = — Finally, 
comparing the other parts of the correlators ([35)1 . (j3T?j) and (j3*Uj) with the final expression 
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for log one can determine the quantities a rs , b rs and u s as functions of "d" which 
remains unfixed. We spare the reader the lengthy details since they are totally technical 
and not very illuminating. We only stress that no further simplification hypothesis is used 
and no expansion is made in neither 1/N nor in any of the coupling constants. In those 
intermediate steps all correlators are treated as if they were exact. 

In order to make the final answer for the bosonic Lagrangian density more familiar we 
have redefined the unfixed parameter d = —Xk 2 /2. In the momentum space we have: 



N 



N 



N 



\r=l 



\s=l 



N 



r ;i - X)k 2 8 aj 3 + Xk 2 gJ + -= e^W A» E B" r 

J v N 



r=l 



N 



+ E 



Xk 2 9 af3 - Xk 2 g. 



aj3 



(43) 



r=l 



Notice that A and A are interpreted as independent gauge fixing parameters such that the 
generating functional (j4*Tj) reproduces the correlation functions (|50~j) . (j35j) and (pffij) with 



C rs = -~7^$rs (44) 

Xk 2 

D = =i 45 

Integration over the vector fields B° will lead to an effective eletromagnetic theory which 
was also studied in j2I] ( see also [21] ) where we have analysed its pole structure and 
concluded that no tachyons appear just like in its two dimensional couterpart f|22[). 

Comparing (|4H|) with the original Lagrangian (1) written in terms of fermionic fields 
we have the bosonization rules for each fermion flavor ( r = 1, 2, • ■ ■ , N ): 



i> T f„<l>r(k) = t^ a k"B^(k), (46) 
-Mk)W + m)Mk) = -fo ^ggj^M ig. (47) 

Similarly to the D = 2 case, we have obtained (|4fij) from the fact that for each fermion 
flavor we have the tensor structure < j^j u > oc a 9^ + bE^ u , which has been derived by 
calculating the fermionic determinant up to the quadratic order. However, this tensor 
structure must hold beyond that approximation ( non-perturbatively ) as a consequence 
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of the Ward identity < j^d u j u >= which must be true also in D = 3 by using a 
gauge invariant regularization. On the other hand, the bosonization rule (|47|) holds as 
a consequence of the quadratic approximation for the fermionic determinant. Taking 
m — > (or z — > — oo ) in the expressions (|27jl we have EEi — »• and n 2 -> 1/(16 V /Z F). 
Plugging back in (J46|) and (|47|) we recover the result of jS] with /3 — 1/4 and 6 1 = in that 
reference which deals with massless fermions. Our rules (|4"H|) and (|47jl also agree with the 
case of massive fermions treated in We stress however, that in both jH] and jT3] only 
free fermions have been considered while our results have been derived for an interacting 
theory with Thirring and eletromagnetic couplings. This indicates that our bosonization 
rules are rather universal at the quadratic level. In the opposite limit of heavy fermions 
(/c 2 /4m 2 ) = z — > in ()26|) we have, at leading order, il 2 — ► and U\ — ► — 1/(4-71") which 
leaves us, see (14 7|) . with a bosonic local and topological Chern-Simons Lagragian for the 
field Bp in agreement with the findings of JU] (see also ^1] and for higher orders ). 

5 Conclusions 

We have derived bosonization maps for the U(l) currents and the fermion Lagrangian 
densities for both QED2 and QED3 with N fermion flavors and Thirring self-interaction. 
Both results hold for finite fermion masses and no derivative expansion on — is made 
as in |3*T] . By turning off the interactions we can reproduce the results of ^3] for a 
particular choice of the regularization parameter in that reference. Our calculations show 
that the U(l) currents, when written in terms of bosons, must be identically conserved 
(topological) as a direct consequence of the Ward-identity < j^d u j u > = which is true 
non-perturbatively if there is no U(l) anomaly. In particular, there is no need for looking 
at higher point functions or compute the fermionic determinant beyond the quadratic 
approximation to confirm that. On the other hand, our maps for the Lagrangian densities 
(j2IJ) and (j4"7J) . which are also independent of the interactions, hold only perturbatively due 
to our quadratic approximation for the fermionic determinant which amounts to neglect 
terms of order higher than (l/iV) 1//2 in the fermionic determinant and consider at most 
two point current correlators to derive the bosonization rules. 

Finally, we should mention that the results derived here could have been obtained 
by a technically more direct way along the lines of [2U| (see also [2H]) but the approach 
used here does not use any auxiliary field and clarifies the fundamental role of the two 
point correlators. Besides, it might be also useful (work in progress) for deriving approx- 
imate bosonic maps for other fermion bilinears like the mass term ipip in D = 3, which 
apparently could not be done by using the approach of |2Uj . 
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